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(57) Abstract: An X-ray imaging method forms a set of a plurality of two-di- 
mensional X-Ray projection images of a medical or veterinary object to be 
examined through a scanning rotation by an X-Ray source viz a viz the object. 
Such X-Ray images are acquired at respective predetermined time instants with 
respect to a functionality process produced by the object. From said set of 
X-Ray projection images by back-projection a three-dimensional volume im- 
age of the object is reconstructed. In particular, an appropriate motion correc- 
tion is derived for the respective two-dimensional images, and subsequently 
as based on a motion vector field from the various corrected two-dimensional 
images the intended three-dimensional volume is reconstructed. 
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Motion-corrected three-dimensional volume imaging method 



BACKGROUND OF THE INVENTION 

The invention relates to an X-ray imaging method., comprising the forming of 
a set of a plurality of two-dimensional X-Ray projection images of a medical or veterinary 
object to be examined through a scanning rotation by an X-Ray source viz a viz said object, 
5 which X-Ray images are acquired at respective predetermined time instants with respect to a 
functionality process produced by said object, and from said set of X-Ray projection images 
reconstructing by back-projection a three-dimensional volume image of said object as being 
recited in the preamble of Claim 1 . 

Such methods have been in use for cardiac and other medical diagnostic 

10 operations. The prime problem for acquiring an appropriate 3D or 4D data set from the 2- 
dimensional images is the movement of items within the human or other body under 
consideration, which movement is sometimes substantially periodic, but not completely so. In 
other circumstances, the object in question, such as an intestine can be considered as moving 
in quite an irregular manner. Furthermore, the object may undergo secondary motion, such as 

1 5 caused by breathing or consciously moving. It is therefore a problem that various two- 
dimensional images donot refer back to an identical 3D object, which would render their 
applicability less good. However, for a top-quality reconstruction the number of processed 
two-dimensional images should be as high as possible, whilst raising the amount of available 
data that is in the same "phase" of the motion, as far as such phase is detectable. 

20 

SUMMARY TO THE INVENTION 

In consequence, amongst other things, it is an object of the present invention 
to raise the number of useful two-dimensional images by providing motion correction for the 
respective two-dimensional images to reconstruct from various corrected two-dimensional 
25 images the intended three-dimensional volume. As disclosed hereinafter in additional detail, 
the correction may be based on feature extraction and/or three-dimensional feature matching, 
such as through markers, bifurcations, or other. 

Now therefore, according to one of its aspects, the invention is characterized 
according to the characterizing part of Claim 1. 
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The invention also relates to a three-dimensional X-Ray apparatus being 
arranged for implementing a method as claimed in Claim 1 . Further advantageous aspects of 
the invention are recited in dependent Claims. 



5 BRIEF DESCRIPTION OF THE DRAWING 

These and further features, aspects and advantages of the invention will be 
discussed more in detail hereinafter with reference to the disclosure of preferred 
embodiments of the invention, and in particular with reference to the appended Figures that 
illustrate: 

10 Figure 1, an X-Ray imaging apparatus, in which the invention may be applied; 

Figure 2, a cardiac cycle and various applicable data acquisition instants 



therein: 



Figure 3, a flow chart of a method embodiment. 



1 5 DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 

Figure 1 illustrates an exemplary X-Ray imaging apparatus in which the 
invention may be applied. Basically, the apparatus allows to form two-dimensional X-Ray 
images of an object to be examined, in particular an object that moves (quasi-)periodically, 
such as a heart and its associated coronary vascular system, or alternatively, an object that 

20 moves rather unpredictably, such as an intestine or part thereof. From the combining of 
various such two-dimensional images, a three-dimensional volume of the object should be 
obtained. By itself, such three-dimensional reconstruction is state of the art. Similar 
technology may be used to produce a four-dimensional data set such as representing a 
selected part of the heart during a phase interval of its motion. 

25 Now, the imaging apparatus 1 includes a C-arm 10 that is mounted on a 

partially shown stand 1 1 . The C-arm can be rotated over an angle such as 180° around its 
center in the direction of double arrow 20 through a motor drive not shown. The C-arm 
accommodates an X-Ray source 12 and an X-Ray image pick-up 13, that are aligned relative 
to each other in such manner that an X-Ray image can be formed of a certain volume around 

30 the above center. These plural of X-Ray images show the volume under examination 

generated by respective different angular orientations of the image forming system 12, 13, 
that are in part shown by dashed lines. The pick-up device may be a series arrangement of an 
X-Ray image intensifier that feeds a television chain, while the signals furthermore are A/D 
converted (14) and stored (15), so that the complete examination yields a series of images 
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(. . .Dm, Di, Di+i ? Di+2 . . .)• These X-Ray images can themselves be processed by known 
reconstruction methods (16) to obtain a three-dimensional volume for examination. This 
volume or various projections therethrough can be displayed on monitor 18. The various 
subsystems of the imaging apparatus are controlled by controller 17. Another applicable 
5 apparatus could be based on multi-slice CT-seanning. 

Figure 2 illustrates a cardiac cycle and various applicable data acquisition 
instants therein. For simplicity, the cardiac cycle H has been shown as ideally periodic; in 
practice, large variations may occur. Preferred data acquisition instants for a three- 
dimensional image are located during end diastole ED on an EGG scale; during this period, 

10 the heart geometry is relatively stationary. The scanning speed of the arm may be 180° in 
some 15 to 20 seconds., at a standard heartbeat rate of some 60 beats per minute, and 
delivering 25 images per second. This would imply an angle variation of some 12-15 degrees 
and a number of images of some 25 for one heartbeat, as has been indicated in the idealized 
Figure. The motion of an intestine or artery would likewise occur on time scales ranging from 

1 5 fractions of a second to several seconds. 

Now, a cardiologist may position a stent delivery catheter provided with 
markers into a region of interest (ROI). By themselves, such markers have been in wide use 
to facilitate interpreting X-Ray images. The markers may show on the X-Ray images as 
relatively darker, or relatively lighter dots with respect to the contrast fluid. Alternatively, no 

20 contrast fluid is used, and the observable contrast is relative to the blood. Next, after 

positioning the center of rotation near the ROI, by moving either the patient or the apparatus, 
a rotational scan of the cardiac ROI is made. Note that the usage of so-called guide wires 
provided with multiple markers allows to use a larger ROI, or rather, to improve the quality 
of the resulting image. 

25 Fundamental in this embodiment is the compensating for cardiac or other 

movement by correcting the two-dimensional images of the rotational scan prior to the back- 
projection. By itself, movement compensation technique has been disclosed in documents EP 
02 292 995.4, EP 03 290 455.9, and EP 02 292 336.1 assigned to the present assignee. Such 
technique allows (1) identifying the markers/features in the two-dimensional image; (2) 

30 motion correction through two-dimensional image-warping through a motion vector field, 
such that identified markers/features are moved to their two-dimensional reference position, 
while using interpolation to correct the rest of the image. Essential is the determination of the 
two-dimensional reference positions for each image to which the markers/features identified 
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in that image are to be corrected to. To determine the two-dimensional reference position of 
each marker/feature in each image, the following steps are executed: 

ECG-based Reference Image selection: the ECG signal is used to select two reference images 
that are in the same cardiac cycle phase and under an appreciable angular distance. Such 
5 angular distance may range between about 90° and for example 45°. Preferably, a relatively 
large rotary angle should be used effectively. Cardiac phase determination via ECG signal 
analysis has been state of the art for both Cardiac Modelling and for Cardiac Reconstruction. 
Reference images with approximately 90° are used in Cardiac Modelling. Other techniques 
for such selection are likewise applicable. 

10 Reference marker/feature Identification in 3D. The markers/features are identified on one of 
the reference images. This identification may need only two markers/features selected by a 
user person, and may even be fully automated. Using the epipolar projection line of the 
markers/features, these are automatically identified in any other reference image. Thereby, 
the three-dimensional positions of the markers/features will be known. 

15 Projection of three-dimensional reference markers/features to 2D. For each of the two- 
dimensional images of the rotational scan, the two-dimensional position of the reference 
markers/features is determined by projecting the position of the three-dimensional 
marker/feature onto the image. For each image, this determines the two-dimensional 
reference positions to which the markers/features found in that image will be corrected. 

20 CORRECTION-BASED PROJECTION IMAGE SELECTION 

The "amount of movement correction" can be measured by for example the 
distance between the identified two-dimensional marker position and the reference two- 
dimensional marker position, or through an ECG analysis, or by a combination of the two 
methods. By selecting images for which the amount of movement correction is less than a 

25 certain value, the image quality aspects of the reconstruction can be influenced: a trade-off in 
the number of images used for reconstruction versus the similarity between images can be 
made. By balancing this value, the image quality of cardiac reconstructions can be 
influenced, and therefore, optimized for its intended (clinical) purpose. 



3 0 THREEDIMENSION AL RECONSTRUCTION 

Via the state of the art Feldkamp Back-Projection Reconstruction, a three- 
dimensional volume is created based on the selected projection images. The volume is then 
visualized for clinical analysis. The image quality of the resulting three-dimensional 
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reconstruction is expected to be clinically usable near the markers/features, and therefore, in 
the ROI. 

OVERLAY OF MULTIPLE THREEDIMENS IONAL RECONSTRUCTIONS 
5 Since "cardiac motion compensation 95 is built into three-dimensional Cardiac 

ROI reconstruction multiple runs of the cardiac region may be made and overlaid 
adequately. This requires the markers to be left at the same position in the patient, while 
making the multiple Cardiac ROI Reconstructions. This allows for example a coronary 
artery, a placed stent, and also a plaque within the ROI to be visualized in a single three- 
1 0 dimensional view. 

FURTHER ASPECTS OF THE INVENTION 

Now in the above, a prime difficulty is the correct determination of the gating window for 
deriving a projection of the moving object. Although in principle such gating window may be 

15 derived from analysing an ECG signal simultaneously with taking the measurements, the 

present inventors have concluded that the analysis for determining the gating instants should 
better be based on analysing the data proper. An advantageous embodiment therefore is 
following one or more feature points, such as arterial bifurcation points or bone elements, 
which would be quite feasible over a heart cycle. Such would then essentially imply a four- 

20 dimensional motion tracking of the coronary feature points, which will be more in detail 
discussed hereinafter. The temporally resolved three-dimensional locations of the feature 
points over the cardiac cycle can be used for velocity and acceleration analysis of the 
coronaries, which would yield optimally suited feature-point-dependent temporal gating 
windows. 

25 One of the particular problems with image reconstruction in a heart 

environment is the limited number of projections that are available for such reconstruction, 
which scarcity will often cause severe streaking artifacts in the final reconstructed image. 
Therefore, the number of usable projections should be substantially increased. This object is 
favourably effected by producing a two-dimensional motion compensation for those 

30 projections that by themselves have not been acquired in the most advantageous cardiac 

phase so that they would be transformed to represent their shape in such most advantageous 
phase. The embodiment approach in particular separates the motion of coronaries into two 
main components: 

a non-linear component caused by the contraction of the heart; 
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a linear component caused by the rotation of the object. 

The result is a transformation matrix that transforms the projection from one 
cardiac phase into another. Now, the motion-compensated two-dimensional projections are 
obtained as follows, whilst referring to Figure 3 that illustrates a flow chart of the method. 
5 Now first, in block 30, the necessary hardware and software facilities are assigned. Then 
subsequently: 

Acquiring a rotational angiography data set from a calibrated system (block 32); 
Reconstructing a low-spatial-resolution volume data set for a specific heart phase (block 34); 
(optional) Estimating a three-dimensional centerline in the volume data (block 36); 
10 Forward projecting the volume data (or rather, the three-dimensional centerline) into the 
successively acquired projections with different projection geometry (block 38); 
Using the forward projected volume, or rather, the centerline (cf. #3) as an initial 
approximation for the correct motion-compensated projection for this viewing angle (block 
40); 

1 5 Calculating a transformation matrix between the initial approximation and the real acquired 
projection at the current viewing angle (block 42); 

Transforming the acquired projection into the correct cardiac phase (block 44); 

Incorporating the additionally acquired projections in the three-dimensional reconstruction 

procedure by the successive application of the above method on any or all appropriate 
20 projections (block 46). Subsequently, the process is stopped, or another series of 

measurements is undertaken (block 48). 

Now, the present invention has hereabove been disclosed with reference to 

preferred embodiments thereof. Persons skilled in the art will recognize that numerous 

modifications and changes may be made thereto without exceeding the scope of the appended 
25 Claims. In consequence, the embodiments should be considered as being illustrative, and no 

restriction should be construed from those embodiments, other than as have been recited in 

the Claims. 
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CLAIMS: 



1. Ail X-ray imaging method comprising the following steps: 

forming a set of a plurality of two-dimensional X-Ray projection images of a 
medical or veterinary object to be examined through a scanning rotation by an X-Ray source 
viz a viz said object, which X-Ray images are acquired at respective predetermined time 
5 instants with respect to a functionality process produced by said object; 

from said set of X-Ray projection images reconstructing by back-projection a 
three-dimensional volume image of said object, 

said method being characterized by deriving an appropriate motion 
correction for the respective two-dimensional images as based on a motion vector field, and 
10 subsequently from the various corrected two-dimensional images reconstructing the intended 
three-dimensional volume. 

2. A method as claimed in Claim 1, wherein said motion correction is derived 
from reference images that are acquired in corresponding instants of the movement of the 

15 object in question that is substantially periodic, and which reference images have 
substantially differing projection orientations. 

3. A method as claimed in Claim 2, wherein said corresponding instants refer to 
corresponding phases of a cardiac movement. 

20 

4. A method as claimed in Claim 3, wherein said movement is derived from 
following one or more feature points of the object, such as bifurcation points. 

5. A method as claimed in Claim 1, and being based on feature extraction for 
25 deriving said motion vector field. 

6. A method as claimed in Claim 1, wherein two-dimensional projections are 
corrected towards their calculated shape in a more advantageous phase of the motion by the 
functionality process of the object. 
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7. A method as claimed in Claim 1, through separating an estimated motion of 
parts of said object into a non- linear temporal component caused by overall contraction 
within said object, and a linear temporal component caused by overall rotation within said 

5 object 

8. A method as claimed in Claim 7, and in particularly applied to coronary 
arteries. 

10 9. A method as claimed in Claim 2, wherein said projection orientations differ by 

an angle in a range between substantially 45 degrees and 90 degrees. 

10. A method as claimed in Claim 1, for use with a coronary artery with a stent in 
place and an artery wall section of said artery being under investigation. 

15 

11. A method as claimed in Claim 1, whilst deriving said motion correction from 
physical elements present in the object, such as markers provided on a stent delivery catheter 
or on a guide wire. 

20 12. A method as claimed in Claim 1, whilst including in said correction an overall 

translation pertaining to said object. 

13. A method as claimed in Claim 1, whilst deriving an amount of movement 
correction from a measured distance between an identified two-dimensional marker/feature 

25 position and a reference two-dimensional marker/feature position, or through an ECG 
analysis, or through a combination of the two methods. 

14. A method as claimed in Claim 1, whilst using built-in cardiac motion 
compensation for three-dimensional cardiac ROI reconstruction, and generating and 

30 overlaying multiple runs of a cardiac region whilst maintaining one or more markers at the 
same position, and by overlaying making the multiple cardiac ROI reconstructions. 

15. A method as claimed in Claim 1, and applied to generating a four-dimensional 
data set. 
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16. A method as claimed in Claim 1, whilst determining a temporal gating as 

being based on a three-dimensional resolving of a feature point location. 

5 17. A method as claimed in Claim 1 , and including one or more steps of the 

following sequence of steps: 

Acquiring a rotational angiography data set from a calibrated system; 
Reconstructing a low-spatial-resolution volume data set for a specific heart phase; 
Estimating a three-dimensional centerline in the volume data; 
10 - Forward projecting the volume data or the three-dimensional centerline into the 
successively acquired projections with different projection geometry; 

- Using the forward projected volume or the centerline as an initial approximation for the 
correct motion-compensated projection for this viewing angle; 

- Calculating a transformation matrix between the initial approximation and the real 
1 5 acquired projection at the current viewing angle; 

Transforming the acquired projection into the correct cardiac phase ; 
Incorporating the additionally acquired projections in the three-dimensional 
reconstruction procedure by the successive application of the above on any or all 
appropriate projections. 

20 

18. An X-Ray apparatus being arranged for implementing a method as claimed in 

Claim 1, and comprising an X-Ray facility for forming a set of a plurality of two-dimensional 
X-Ray projection images of an object to be examined through a scanning rotation by an X- 
Ray source vis a vis said object, which X-Ray images are acquired at respective 
25 predetermined time instants with respect to a functionality process produced by said object; 

data processing means fed by said X-Ray facility for from said set of X-Ray 
projection images reconstructing by back-projection a three-dimensional volume image of 
said object, 

and correcting means interacting with said data processing means for by 
30 deriving an appropriate motion correction for the respective two-dimensional images as 
based on a motion vector field, and for subsequently feeding the various corrected two- 
dimensional images to said data processing means for reconstructing the intended three- 
dimensional volume. 
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